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We calculate the cross section of inclusive direct J/ip photoproduction at next-to- leading order 
within the factorization formalism of nonrelativistic quantum chromodynamics, for the first time 
including the full relativistic corrections due to the intermediate ^s}^\ ''^S^^\ and '^fj*' color-octet 
states. A comparison of our results to recent HI data suggests that the color octet mechanism is 
indeed realized in J/tp photoproduction, although the predictivity of our results still suffers from 
uncertainties in the color-octet long-distance matrix elements. 
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The factorization formalism of nonrelativistic quan- 
tum chromodynamics (NRQCD) [l| provides a consis- 
tent theoretical framework for the description of heavy- 
quarkonium production and decay, which is known to 
hold through two loops 2]. This implies a separation 
of process-dependent short-distance coefficients, to be 
calculated perturbatively as expansions in the strong- 
coupling constant a^, from supposedly universal long- 
distance matrix elements (LDMEs), to be extracted from 
experiment. The relative importance of the latter can 
be estimated by means of velocity scaling rules; i.e., the 
LDMEs are predicted to scale with a definite power of the 
heavy-quark (Q) velocity v in the limit v ^ 1. In this 
way, the theoretical predictions are organized as double 
expansions in as and v. A crucial feature of this formal- 
ism is that it takes into account the complete structure 
of the QQ Fock space, which is spanned by the states 
n = 2S+i^M ^j^j^ definite spin S, orbital angular mo- 
mentum L, total angular momentum J, and color multi- 
plicity a = 1, 8. In particular, this formalism predicts the 
existence of color-octet (CO) processes in nature. This 
means that QQ pairs are produced at short distances in 
CO states and subsequently evolve into physical, color- 
singlet (CS) quarkonia by the nonperturbative emission 
of soft gluons. In the limit v ^ 0, the traditional CS 
model (CSM) is recovered in the case of S-wave quarko- 
nia. In the case of J/tp production, the CSM prediction 

is based just on the CS state, while the leading rel- 
ativistic corrections, of relative order 0{v'^), are built up 
by the ^4^', ^sf \ and ^pf (J = 0, 1, 2) CO states. 

Fifteen years after the introduction of the NRQCD fac- 
torization formalism the existence of CO processes 
and the universality of the LDMEs are still at issue and 
far from proven, despite an impressive series of exper- 
imental and theoretical endeavors. The greatest suc- 
cess of NRQCD was that it was able to explain the J/ip 
hadroproduction yield at the Fermilab Tevatron Q , while 
the CSM prediction lies orders of magnitudes below the 
data, even if the latter is evaluated at next-to-leading 
order (NLO) or beyond [1, 0| . Also in the case of J/ip 
photoproduction at DESY HERA, the CSM cross section 



significantly falls short of the data, as demonstrated by 
a recent NLO analysis [H using up-to-date input param- 
eters and standard scale choices, leaving room for CO 
contributions 0. Similarly, the J/ip yields measured in 
electroproduction at HERA and in two-photon collisions 
at CERN LEP2 were shown [s','^ to favor the presence of 
CO processes. As for J/tp polarization in hadroproduc- 
tion, neither the leading-order (LO) NRQCD prediction 
[lol |. nor the NLO CSM one 0] leads to an adequate 
description of the Tevaton data. The situation is quite 
similar for the polarization in photoproduction at HERA 
i. 

In order to convincingly establish the CO mechanism 
and the LDME universality, it is an urgent task to com- 
plete the NLO description of J/^jJ hadro- [1, [^, [ll|l and 
photoproduction II, [ill, regarding both J/tp yield \4.\l^ 
and polarization [5|, |6|, |ll[ , by including the full CO con- 
tributions at NLO. While the NLO contributions due to 
the ^sl^^ and ^S^i^ CO states may be obtained [11] using 
standard techniques, familiar from earliest NLO CSM 

calculations [l^, the NLO treatment of ^Pf^ states in 
2—^2 processes requires a more advanced technology, 
which has been lacking so far. In fact, the "^pj^' contri- 
butions represent the missing links in all those previous 
NLO analyses [^-iS,, Jd,, ,12], and there is no reason at 
all to expect them to be insignificant. Specifically, their 
calculation is far more intricate because the application 
of the ^Pf^ projection operators to the short-distance 
scattering amplitudes produce particularly lengthy ex- 
pressions involving complicated tensor loop integrals and 
exhibiting an entangled pattern of infrared (IR) singu- 
larities. This technical bottleneck, which has prevented 
essential progress in the global test of NRQCD factor- 
ization for the past fifteen years, is overcome here for 
the first time. So far, only two complete NLO analyses 
of heavy-quarkonium production in high-energy collisions 
involving CO states have beenperformed: the total cross 
section of hadroproduction 13] and the inclusive cross 
section at finite transverse momentum px in two-photon 
collisions [l3|. However, the former case corresponds to 
a 2 — > 1 process, which enormously simplifies the calcu- 
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FIG. 1: Sample diagrams contributing at LO (a and d) and 
to the virtual (b and e) and real (c and f) NLO corrections. 
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FIG. 2: Overview of the IR singularity structure. 

lation, and the latter case does not involve virtual cor- 
rections in P-wave channels. 

In direct photoproduction, a quasi-real photon 7 that 
is radiated off the incoming electron e interacts with a 
parton i stemming from the incoming proton p. Invoking 
the Weizsacker- Williams approximation and the factor- 
ization theorems of the QCD parton model and NRQCD 
[H, the inclusive J/-0 photoproduction cross section is 
evaluated from 

da{ep ^ J/i) + X) = dxdy f^/e{x)f,/p{y) (1) 

X {0''^'^[n])d(j{ji cc[n] + X), 

where f-^/eix) is the photon flux function, fi/p{y) are 
the parton distribution functions (PDFs) of the proton, 
{0-^^'f'[n]) are the LDMEs, and da{ji cc[n] + X) are 
the partonic cross sections. Working in the fixed-flavor- 
number scheme, i runs over the gluon g and the light 
quarks q = u,d, s and anti-quarks q. 

The Feynman amplitudes of 71 — >■ cc[n\ + X are cal- 
culated by the application of appropriate spin and color 
projectors onto the usual Feynman amplitudes for open 
cc production (isj . Example Feynman diagrams for 
partonic LO subprocesses as well as virtual- and real- 
correction diagrams are shown in Fig.[T] Important prop- 
erties of these projections are that the relative momen- 
tum q between the c and c quarks has to be set to zero, 
in the case of P-wave states after taking the derivative 
with respect to q. 

We checked analytically that all appearing singulari- 
ties cancel. As for the ultraviolet singularities, we renor- 
malize the charm-quark mass and the wave functions of 
the external particles according to the on-shell scheme, 
and the strong-coupling constant according to the mod- 
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FIG. 3: Separate and joint dependences of da{ep J/^ + 
X)/dp^ at = 20 GeV^ in full NRQCD at LO and NLO on 
and ^f. 



ificd minimal-subtraction scheme. Figure [5] displays an 
overview of the IR singularity structure. In the case of 
the ^S\^\ ^S^\ and '^sf^ states, the soft and coUinear 
singularities of the real corrections are canceled as usual 
by complementary contributions stemming from the vir- 
tual corrections and by the absorption of universal parts 
into the proton and photon PDFs, the latter entering via 

resolved photoproduction. In case of the ^P^^ states, the 
soft singularity structure is more complex. The reason is 
the following: In the soft limit, the real-correction ampli- 
tudes factorize into LO amplitudes and so-called eikonal 
factors. Taking the derivative with respect to q and 
squaring the amplitudes then leads to additional soft #2 
and soft #3 terms because the derivative has to be taken 
of the eikonal factors as well. The soft #3 terms are pro- 
portional to a linear combination of the short-distance 
cross sections to produce the and 35^1 states. They 
are canceled against IR singularities stemming from ra- 
diative corrections to the (©'^/'^(^S'i^l)) and (©"'/'^(^S'fi)) 
LDMEs. The soft #2 terms do not factorize to LO cross 
sections. They also cancel against virtual-correction con- 
tributions as the usual soft jfl terms. 

Apart from the analytical cancellation of all occur- 
ring singularities, our calculation passes a number of fur- 
ther nontrivial checks. We implemented two indepen- 
dent methods for the reduction of the tensor loop inte- 
grals, which yielded identical results. As for the real cor- 
rections, the numerical evaluation of our expressions for 
the squared matrix elements agree with numerical out- 
put generated using the program package MadOnia [l^ , 
well within the numerical uncertainty of the latter. We 
verified that our results are stable w.r.t. variations of the 
phase space slicing parameters introduced as a demarca- 
tion between the soft and/or coUinear regions from the 
rest of the three-particle phase space. [? ]. Finally, we 
could nicely reproduce the NLO CSM results of Ref. [11] 
after adopting the inputs chosen therein. For space lim- 
itation, we refrain from presenting here more technical 
details, but refer the interested reader to a forthcoming 
publication. 
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We now describe our theoretical input and the kine- 
matic conditions for our numerical analysis. We set 
™c = adopt the values of raj/^, rrie, and 

a from Ref. [16|, and use the one- loop (two- loop) for- 
mula for ai"^^(/x), with = 3 active quark flavors, 
at LO (NLO). As for the proton PDFs, we use set 
CTEQ6L1 (CTEQ6M) [17] at LO (NLO), which comes 

with an asymptotic scale parameter of Aq^j-, = 215 MeV 

(326 MeV), so that A^^^^^ = 249 MeV (389 MeV). 
We evaluate the photon flux function using Eq. (5) of 
Ref. fTsi with the cut-off Ql^^^ = 2 GeV^ [11 [3 on the 
photon virtuality. Our default choices for the renormal- 
ization, factorization, and NRQCD scales are fi^ — l^f ^ 
rriT and /iA = Wc, respectively, where rriT = y^|r+4^7if 
is the J /ij] transverse mass. We adopt the LDMEs from 
Ref. 21], which were fitted to Tevatron I data using 
the CTEQ4 PDFs, because, besides the usual LO set, 
they also comprise a higher-order-improved set deter- 
mined by approximately taking into account dominant 
higher-order effects due to multiple-gluon radiation in in- 
clusive J /if) hadroproduction, which had been found to 
be substantial by a Monte Carlo study 22]. This obser- 
vation is in line with the sizable NLO corrections recently 
found in Refs. 

0,11 [HI- still excluding the ^Pf^ channels 
at NLO. Of course, LDME fits to more recent Tevatron 
data are available, but their goodness is clearly limited by 
the present theoretical uncertainties in the short-distance 
cross sections, preventing the increase in experimental 
precision gained since the analysis of Ref. [211] from actu- 
ally being beneficial. Apart from that, the central values 
of the J/tp LDMEs have only moderately changed, as 
may be seen by comparing the LO results of Ref. [2l| 
with those recently obtained [l^ by fitting Tevatron II 
data using the CTEQ6L1 PDFs Because the pr 

distributions of the ^S'q*' and ^Pj*' contributions to J/ij} 
hadroproduction exhibit very similar shapes, fits usually 
only constrain the linear combination 

with an r value of about 3.5 [2l|, [2^. As in 
Ref. Q, we take the democratic choice {0-^''l'(^sf^)) = 
(r/m2)(0'^/'/'(3Ff|*l)) = as our default. 

Recently, the HI Collaboration presented preliminary 
data on inclusive J /ij) photoproduction taken in colli- 
sions of 27.6 GeV electrons or positrons on 920 GeV pro- 
tons in the HERA II laboratory frame [l^]. They ni cely 
agree with their previous measurement at HERA I [19|. 
These data come as singly differential cross sections in 
Pt, W = V(P7 +PpY^ and z = {pjj^ ■Pp)/{p-i ■ Pp), in 
each case with certain acceptance cuts on the other two 
variables. Here, Pj, Pp, and pj/^ are the photon, proton, 
and J/'tp four-momenta, respectively. In the comparisons 
below, we impose the same kinematic conditions on our 
theoretical predictions. 

We start our numerical analysis by estimating the theo- 
retical uncertainties. The dependences on the unphysical 



scales and ^if are investigated in full NRQCD at LO 
and NLO for the typical case of da{ep ^ J/ip -\- X)/dp^ 
at p^ = 20 GeV^ in Fig. |3l Contrary to naive expec- 
tations, the scale dependence is not reduced when pass- 
ing from LO to NLO. Detailed investigation reveals that 
this behavior may be ascribed to the fact that the new 
coefficient of al{^r) is greatly dominated by the part 
that does not carry logarithmic dependence on fir or 
/i / , mainly arising from the gluon- induced ^s]^^ and ^Pq^^ 
channels, while the complementary part still formally 
warrants renormalization group invariance up to terms 
beyond NLO. As for the dependence on mc, a reduction 
of rric from mj/^/2 « 1.55 GeV to 1.4 GeV typically 
entails a rise in cross section by about 50%. The free- 
dom in sharing m/^"^ of Eq. © between (©"'/'^(i^*')) 
and (r/TO^)(0"'/^("^PQ^')) typically creates an uncertainty 
of about 10%. The bulk of the theoretical uncertainty 
is actually due to the lack of knowledge of the com- 
plete NLO corrections to the cross section of inclusive 
J/ip hadroproduction, which is instrumental for a reli- 
able NLO fit of the CO LDMEs to the Tevatron data. 
As explained above, these corrections are expected to be 
dominated by positive and sizable contributions from real 
QCD bremsstrahlung [1, |^, [ll|, [ll], leading to a signifi- 
cant reduction of the CO LDMEs [2l|. At present, the 
theoretical uncertainty in inclusive J/ijj photoproduction 
from this source may be conservatively estimated by com- 
paring the full NRQCD evaluations using the LO and 
higher-order- improved LDME sets of Ref. [2l| , with the 
understanding that the former is bound to overshoot a 
future evaluation with a genuine NLO set. This kind 
of uncertainty is indicated in the remaining figures by 
shaded (yellow) bands, whose upper margins (solid lines) 
refer to the LO set. 

The HI measurements [l^ [2^ of the p^, W, and z 
distributions of inclusive J ftp photoproduction are com- 
pared with our new NLO predictions in full NRQCD in 
Fig. |31[a)-(c), respectively. For comparison, also the de- 
fault predictions at LO (dashed lines) as well as those 
of the GSM at NLO (dot-dashed lines) and LO (dotted 
lines) are shown. Notice that the experimental data are 
contaminated by the feed-down from heavier charmonia, 
mainly due to ijj' — J/ip -\- X , which yields an estimated 
enhancement by about 15% jl2| . Furthermore, our pre- 
dictions do not include resolved photoproduction, which 
contributes appreciably only at z < 0.3 [2lj, and diffrac- 
tive production, which is confined to the quasi-elastic do- 
main at z w 1 and px ~ 0. These contributions are effi- 
ciently suppressed by the cut 0.3 < z < 0.9 in Figs.llja) 
and (b), so that our comparisons are indeed meaning- 
ful. We observe that the NLO corrections enhance the 
NRQCD cross section, by up to 115%, in the kinematic 
range considered, except for z < 0.45, where they are 
negative. As may be seen from Fig. Hl^c), the familiar 
growth of the LO NRQCD prediction in the upper end- 
point region, leading to a breakdown at z = 1, is further 
enhanced at NLO. The solution to this problem clearly 
lies beyond the fixed-order treatment and may be found 
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FIG. 4: (a) (b) W , and (c) z distributions of inclusive J/^ photoproduction at LO and NLO in the GSM and full NRQGD 
in comparison with HI data [l9l. [201]. The shaded (yellow) bands indicate the theoretical uncertainty due to the GO LDMEs. 



in soft collinear effective theory [2j|. Ttie experimental 
data are nicely gathered in the central region of the er- 
ror bands, except for the two low-z points in Fig. |4jc), 
which overshoot the NLO NRQCD prediction. However, 
this apparent disagreement is expected to fade away once 
the NLO-corrected NRQCD contribution due to resolved 
photoproduction is included. In fact, the above consid- 
erations concerning the large size of the NLO corrections 
to hadroproduction directly carry over to resolved pho- 
toproduction, which proceeds through the same partonic 
subprocesses. On the other hand, the default CSM pre- 
dictions significantly undershoot the experimental data, 
by typically a factor of 4, which has already been ob- 
served in Ref. Q. Except for > 4 GeV^, the situation 
is even deteriorated by the inclusion of the NLO correc- 
tions. 

Despite the caveat concerning our limited knowledge 
of the CO LDMEs at NLO, we conclude that the HI 



data [19|, |20| show clear evidence of the existence of CO 
processes in nature, as predicted by NRQCD, supporting 
the conclusions previously reached for hadroproduction 
at the Tevatron Q and two-photon collisions at LEP2 
0. In order to further substantiate this argument, it is 
indispensable to complete the NLO analysis of inclusive 
J/V' hadroproduction in NRQCD, by treating also the 

■^pj*' channels at NLO, so as to permit a genuine NLO fit 
of the relevant CO LDMEs to Tevatron and CERN LHC 
data. This goal is greatly facilitated by the technical 
advancement achieved in the present analysis. 
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